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I.    INTRODUCTION 

A.   BACKGROUND 

The  art  and  science  of  marine  positioning  date  back  to  antiquity, 
and   have  progressed  from  black  magic,    astronomical   positioning,    radar, 
Loran   C  and  Omega  to  artificial   satellite  systems    (Bresslau,    1980,    pp. 
A35) .      The  accuracy  of   land-based   electronic  systems   varies:      e.g., 
the   Decca    Pulse-8  system   has   an   accuracy  of   ±   50  m   at  200   nmi; 
Loran-C   has  an  accuracy  of  ±  900  m  at  2000   km   (Laurila,    1976,    pp. 
348,429);    Omega   has   an   accuracy  of  ±   7.2   km  globally    (Creamer  et  al., 
1985, pp.    269);    while  the   Navy   Navigational    Satellite  System   has   a  30-m 
accuracy  with   precise  ephemerides   for  a   single  pass.      In   general   the 
accuracies  of  electronic  positioning   systems   are  degraded   by   land 
propagation   distortion    (Laurila,    1976,    pp.    502). 

Astronomical    navigation   methods   have  serious  observational    limita- 
tions,   i.e,    the  measurement  depends   on   a   plumb   line  and   it   is   weather- 
sensitive,    although   a   global   datum   is   possible.      The  concept  of 
positioning   by  an   artificial    satellite   system   began   with   the   use  of  the 
Baker-Nunn   camera  for  direction   determination.      This   formed  the  foun- 
dation  of   real   3-dimensional   positioning    (Laurila,    1976,    pp.    383). 

The  current  operational    navigation   satellites   are  part  of  the   U.S 
Navy   Navigation    Satellite  System,    usually  called   TRANSIT.      The  system 
was   conceived   in    1958  at  the  Applied    Physics    Laboratory  of  the  Johns 
Hopkins   University  and   became  operational    in  January   1964  and   commer- 
cially  available  in    1967.      The  technical   approach   was   to   use  the   Doppler 
shift  to  provide   intermittent  position   fix    updates    (Laurila,    1976,    pp. 
463).      The  TRANSIT   system  consists  of  five  operational   satellites   in 
polar  orbits  with   altitudes  of   1100   km,    periods   of   105  min,    and   frequen- 
cies of   150  MHz  and  400  MHz.      The  satellites   continuously   broadcast 
predicted  orbital   coordinates   along   with   timing   signals,    and   are  checked 
by  four  monitoring   stations   and  two   injection   stations    (Laurila,    1976, 
pp.    464). 


In   addition,    the   Department  of   Defense   has   been   developing  the 
Navigation   Satellite  Timing   and    Ranging  Global    Positioning   System 
(NAVSTAR  GPS)    under  the  Joint   Program  Office   (JPO)    since   1973  as   a 
merger  of  Navy's  Time  and   Navigation    (TIMATION),    with   Air   Force 
Project  621 B    (Boissenin,    1985,    pp.    1).      The  basic  method  of 
GPS/TRANSIT  position   fixing   differs  from  conventional   celestial   naviga- 
tion or  visual   satellite  observation   in   that   range,    range-rate  and  angle 
are  measured,    since  electronic  systems   do  better  in    ranging  than   in 
angle  measurement   (Mueller,    1964,    pp.    235-268;    Easton,    1980,    pp.    14). 
The  TIMATION    is  a   program  to  advance  the  development  of   high 
stability  oscillators,    time  transfer,    and  two-dimensional    navigation; 
Project  621b   is  to  develop  a   highly  accurate  3-dimensional    navigation 
system   (Boissenin,    1985,    pp.    1). 

The  development  approach   for  GPS   selected   by   the   Defense   System 
Acquisition    Review  Council    (DSARC)    was   a   step-wise  development  and 
test  program,    leading   in   successive  phases  to  an   operational   GPS.      Each 
phase  was   designed   to  build   and   expand  on   the  previous   phase   in   an 
integrated   and   cohesive  manner.      The  decision   of   DSARC    I    in    December 
1973  was   to  proceed  with   a   phase   I    Concept  and  Validation   to  concen- 
trate on    validation   of  design   concepts   through    Development  Test  and 
Evaluation    (DT&E)    of   user  equipment.      The  decision   of   DSARC    11    in 
June   1979   approved   proceeding    into  full    scale  engineering   development 
to  complete  the   DT&E  and    Initial   Operational   Test  and   Evaluation 
(lOT&E)   of   user  equipment.      A   DSARC    III    decision    is   planned  for 
October   1985  to  obtain    Approval   for   Full    Production    (AFP)    of  the   user 
equipment    (Boissenin,    1985,    pp.    2). 

B.      GPS   POSITION   AND   NAVIGATION   TECHNIQUE 

The  Global    Positioning   System    (GPS)    is   a   space-based    radio   posi- 
tioning  navigation   system  that   is   designed   to   provide  continuous,    highly 
accurate,    three-dimensional   position   data   to  within   a    16-m   spherical 
error  probability    (SEP),    velocity  to  within   0.05-0.10  m/s   and   system 
time  to  within   55   ns  to  suitably  equipped   users   anywhere  within   500   nmi 


10 


of  the  earth.      The  position   solution   is  based  on  World  Geodetic  System 
(WGS)    1972  ellipsoids   in   the  Earth-Centered,    Earth-Fixed    (ECEF)    coor- 
dinates.     Accurate  position   and   velocity  determination   will   permit  GPS  to 
supplant  less  accurate  systems   such   as   Loran   C,    Omega,    and   Decca 
Pulse-8.      GPS   user  equipment   receives   L-band   signals  from  four  satel- 
lites to  measure  pseudoranges   and    range-rates   to  solve  for  positions    in 
three  dimensions,    velocity  and  time   (Boissenin,    1985,    pp.    2-3). 

C.      GPS  CONFIGURATION 

The  GPS   consists  of   space  segments,    control    segments   and   user 
segments: 

1 .      Space  Segment 

The  space  segment  is   composed  of  a   constellation   of   satellites 
that  continuously  transmit  signals   to  all    users.      When   fully  operational 
the  GPS   satellite  constellation   will   deploy   18  satellites    (plus  3   spares)    in 
circular  orbits   at  altitudes  of   10900   nmi,    inclined  45-65  deg,    and   having 
approximately   12-hr  periods.      The   spacing   of  the  satellites    in   their 
orbits  will   be  arranged   such   that  a   minimum  of  four  satellites   will    be 
available  to  a    user  anywhere  on   the  globe   (Boissenin,    1985,    pp.    2; 
Jorgensen,    1984,    pp.    1).      Operational   GPS   satellites   have  projected 
mean   mission   durations  of  six   years   and   design    lives  of  seven   and   a 
half  years    (Senus   and   Heuerman,    1983,    pp.    73). 

The   unperturbed  orbit  of  satellites   may   be  given   by   six 
Keplerian   elements,    namely  the  semi-major  axis    (a),    eccentricity    (e), 
mean   anomaly    (M)    or  time  of  perigee  passage   (T),    right  ascension   of 
ascending   node   (fl),    inclination   of  orbital   plane  with    respect  to  the 
equatorial   plane   (I),    and   argument  of   perigee   (w) .      The  first  three 
parameters   determine  the  position   of  the  satellite   in    its  orbital   plane, 
and  the   last  three  parameters   determine  the  position   of  the  orbital    plane 
in  the  celestial   coordinate   system    (Torge,    1980,    pp.    115-119;    Mueller, 
1964,    pp.    147-157). 
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The  actual   orbit  of  a   satellite  departs  from  the   Keplerian 
model   due  to  the  effect  of  various   disturbing  gravitational   forces   and 
non-gravitational   forces.      The  orbit  can   be  viewed   as  oscillating   ellipses 
which   are  given   at  each   instant  by  the  current  orbital   elements.      The 
main   gravitational   forces   are  due  to  the  earth,    moon   and   sun.      The 
earth's  gravitational   force  depends  on   satellite  position   and   time  due  to 
irregularities   in   the  earth's   density.      The   non-gravitational   forces   are 
due  to  solar   radiation   pressure,    magnetic  effects,    atmospheric  drag, 
and   relativity.      The  actual   orbit  can   be  modeled   by  the   initial    six 
Keplerian   elements   plus  time   (t),    ellipsoidal    harmonic  coefficients    (zonal 
harmonic  J,    tesseral    harmonics  J   and    K),    and  their   rates  of  change 
(Torge,    1980,    pp.    116-117;    Heiskanen   and   Moritz,    1967,    pp.    342-345). 
The  computation   of  satellite  position   by   ellipsoidal    harmonic  coefficients 
needs   sophisticated   software   in   the   receiver   (MacDoran   et  al.,    1984, 
pp.    65). 

The  actual   orbit  can   be  modeled   also  by  the   initial   six 
Keplerian   elements,    secular  drift  terms   and   harmonic  coefficients   to 
reduce  computation.      Such   a   model,    however,    requires   continuous 
updating.      As   applied    in   the  current  GPS,    corrections   are   uploaded   to 
satellites   every  24   hours   by  the  master  station,    and  these  are   used   to 
make   hourly  changes  on   board  the  satellite   (MacDoran   et  al.,    1984,    pp. 
65). 

GPS   satellites   simultaneously  transmit   navigation    information   on 
two   radio  frequencies,    LI    and    L2.      The   LI    signal    is   modulated   with    C/A 
code    (C/A   =   coarse-acquisition)    and   a    P-code    (P   =   precission),    and   the 
L2   signal    is   modulated  with   the   P-code  only.      The  functions   of  the 
codes   are  to   identify  the  particular   satellite  and   to  measure   phase 
shifts;    both   are   necessary  to  measure  pseudoranges .      Both   the   LI    and 
L2   signals   are  also  continuously   modulated   with   the   navigation   data   bit 
stream  at  50  bps    (Remondi,    1984,    pp.    2;    Milliken   and   Zoller,    1980,    pp. 
6).      Due  to  different  gravitational    potentials   on   the  ground   and   at  the 
satellite,    and   to  velocities  of   users    relative  to  the  satellites,    the  satel- 
lites  actually  transmit  a   slightly   lower  frequency    (Af=  4.45   x    10"^°    so 
that   10.22999999545  MHz   is   shifted  to  a   nominal    10.23  MHz    (Spilker, 
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1980,    pp.   33).      The  navigation  message  contains   information  on  the 
status  of  the  satellite,   the  time  synchronization   information  for  transfer 
from  C/A-  to  P-code,   the  parameter  for  computation  of  the  clock 
correction,   the  ephemerides  for  computation  of  satellite  position,   the 
parameter  for  propagation  delay  correction,    and  the  almanac  for 
computing  Geometric   Dilution   of   Precession    (GDOP)/sateliite  selection. 
The  navigation  message  is   intended  to  keep  the  satellite  ephemerides 
and  timihg   system  as   accurate  as   possible    (Milliken   and   Zoller,    1980, 
pp.    7;   Van   Dierendonck  et  al.,    1980,    pp.    576).      The  contribution  of 
satellite  segment  errors  was   monitored   during  the  GPS   development 
phase,    and  the   recorded   bias   errors  were  found  to   be  for  X,    -2  to   -4 
m;    for  Y,    -5  m;    and  for  Z,    ±3  m    (JPO,    1984,    pp.    6). 

2.  Control  Segment 

The  control   segment  tracks   all   satellites    in   orbit.      The 
purpose  of  satellite  tracking    is   to  observe  the  position   and   motion   of  a 
satellite  as   a   function   of  time  with   a   precision   sufficient  to   refine   its 
equation   of  motion   to  permit   predicting  future   position   with   at   least  the 
accuracy  achievable  by  direct  observation.      The  operational   control 
segment  consists   a   master  control    station,    a   ground   control    station, 
monitor  stations   and   an   alternate  control   station.      The  monitor  stations 
collect  atmospheric  and   satellite  tracking   data  from  all    satellites   in    view, 
which   are  encrypted   and   transmitted   to  the  master  control   station,    the 
ground  control   station   and  the  alternate  control   station    (Boissenin, 
1985,    pp.    2).      The  master  control    station    uses   one  week's   data   from   ail 
monitor  stations  to  control   the  satellites'    health,    to  generate   new   satel- 
lite ephemerides   and   clock  offset   updates,    and   then   transmit  these 
updates   to  the  constellation   of  satellites   every  24   hours    (Milliken   and 
Zoller,    1980,    pp.    69). 

3.  User  Segment 

The  fundamental   capabilities   of  the  GPS   user  segment  are 
satellite  selection,    signal   acquisition,    tracking   and   measurement,    data 
recovery,    and   data   processing   for  the   real   time   user.      The   hardware 
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component  is   comprised  of  an   antenna,    a    receiver,    and   input/output 
devices    (Glazer,    1980,    pp.    81). 

If  a   GPS   receiver  has  a   very  accurate  clock  synchronized 
precisely  with   the  satellite  clock,    it  can   measure  the  transit  time  of  the 
signal   between   the  satellite  and   the   receiver  and  multiply  this   by  the 
speed  of   light  to  get  the   range  from   receiver  to  satellite.      Due  to 
expense  and   the   impracticality  of  maintaining  a   synchronized   clock  at 
both   the  satellite  and  the   receiver,    only  a   crystal   clock   is    used   in   the 
receiver,    so  that  the  measured   signal   transit  time   is   subject  to   receiver 
clock  error.      This   uncorrected   distance   is   called  the  pseudorange.      By 
measuring  four  pseudoranges  to  different  satellites,    the   user's   position 
can   be  calculated   in   the   ECEF  World  Geodetic  System   1972,    which   may 
then   be  transformed   to  another  datum    (Milliken   and   Zoller,    1980,    pp. 
3-4). 

As   the  satellite  signals   vary,    one  can   exploit  each   signal   by 
various   means   to  get  a   position.      At  present  there  are  two  methods: 
The  first   uses   the  code  signal   to  measure  pseudorange  and   pseudorange 
rate;    and   the  second   uses   the   uncoded  transmission   signal   as   a   source 
for  very-long-base-line  interferometry    (VLBI).      Boissenin    (1985), 
Buennagel   et  al.,    (1984),    Bock  et  al.,    (1984),    Borel    (1980),    MacDoran 
(1984),    JPO    (1984),    and  Glazer   (1980)    provide  detailed   information   on 
the   uncoded   GPS    receivers.    . 

The  Texas    Instruments   TI-4100  GPS    receiver  uses   both 
C/A-code  and    P-code  systems.      It   is   a   single-channel,    two-frequency, 
digitized   multiplexing    receiver,    capable  of  tracking   simultaneously   up  to 
four   satellites   with   continuous    LI    and    L2   frequency   measurement. 
Lock-on   for  each    signal    is   for  5  ms,    for  a   total   of  40  ms   to  complete 
one  sequence    (Johnson   et  al.,    1984,    pp.    61).      The   receiver  clock 
permits   connection   to  an   external   atomic  clock   so  that  the  minimum 
number  of  satellites    needed   is    reduced   to  three.      The  standard    receiver 
unit  consists   of  an   antenna   and   pre-amplifier,    a   control   display   unit 
and   a   dual   cassette   recorder.      The  antenna   assembly   consists   of  a 
conical   omni-directional   antenna   and  the  pre-amplifier.      The  Control 
Display   Unit   is   a    hand-held    input/output  device  with   alphanumeric 
displays    (TI-4100,    1983,    pp.    1-4  to   1-7). 
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Data  acquired  with  a  TI-4100  receiver  were  used   in   this   study 
which   involved   processing   pseudorange  data   into  position   data   and 
comparing   GPS-determined   positions  with   Trisponder-determined   position 
data.      The  field  data  were  acquired  by  the  U.S   Naval  Oceanographic 
Office   (NAVOCEANO)    in   Mississippi   Sound.      The  development  of  the 
computer  programs   was   very  slow  due  to  a   large  data   set  that  was   very 
difficult  to  handle;    however,    a   mathematical   model  was   developed. 
Calculations  were  actually  made  at  the   Naval   Surface  Weapon    Center 
(NSWC)    using   NSWC  computer  programs   written   by   S.L.    Meyerhoff 
because  the   NSWC  computer  programs  were  found  to   be   incompatible 
with   the   NPS    IBM-3033;    nor  was    it  possible  to  get  the   NSWC   computer 
program  to   run   on   the   Fleet   Numerical   Oceanography   Center's 
Cyber-175.      The   NSWC   least   squares   adjustment  program  was   more 
sophisticated   than   the  one   I    developed   in   that   it   included   a   combination 
of  bias   Doppler  data,    bias   frequency  of  the  satellite's   clock,    and   bias 
frequency  of  the   receiver's   clock. 
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II.    MATHEMATICAL  MODEL   USED    IN   THE   POINT   POSITIONING   MODE 

The  mathematical   model   may   be  divided   into  two  parts,    namely 
determination   of  coordinates   in  WGS   1972,    and  transformation   to  geodetic 
coordinates   in   some  other  datum.      For  this   evaluation,    the  comparison 
between    launch   positions   determined   by  Trisponders   and   by   GPS   is 
done   in    UTM,    so  a   transformation   of  coordinates   is   performed   for  the 
shore  station. 

To  get  good   accuracy  the   relative  geometry  of  the  satellites   and 
the  ground    receiver  must  be  good.      This   geometry   is   calculated   by  the 
Geometric   Dilution   of   Position    (GDOP)    method.      The   unknown   position   of 
the   receiver's   antenna    (X,    Y,    Z   coordinates)    and    receiver's   clock  offset 
are  functions  of  four  satellite  positions   and  four  corrected   pseudo- 
ranges.      The  satellites'   positions   are  calculated   from  ephemerides.      The 
receiver  measures  four  observed   pseudoranges  to  different  satellites 
which   must  be  corrected   for  atmospheric  delay,    relativity,    earth    rota- 
tion  and   satellite  clock  offset.      The  four   unknown   parameters  of  the 
receiver  are  calculated   either   by   the  conventional    iterative   least   squares 
method  or  the   Kalman   filter  method.      Meyerhoff   (1985),    has   developed 
computer  programs   to  calculate  either   by   least   squares   adjustment  or 
Kalman   filter  methods.      The  formulas    used   for  the  computations    in   the 
computer  program  will   be  mentioned  only   here;    their  derivations   can   be 
found   in   the   references. 

A.      SATELLITE   SELECTION 

If  more  than   four   satellites   are  above  the   user's    horizon,    the 
receiver  can   select  the  best  geometry  to  get  the  best   position,    other- 
wise  satellite   selection    is   made  on   the   basis   of  the  satellites'    health   and 
the  ages  of  the  data   messages.      The  geometric  dilution   of   position 
(GDOP)    is   a   measurement  of  the   best  geometry,    which    is   a   function   of 
the   user's   and  the  satellite's   position.      The  satellite's   coordinates   are 
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computed  from  the  almanac.      The  parameters  of  the  almanac  are  simpli- 
fied  ephemerides   parameters   and   are  stored   in   the   receiver's   memory  or 
can   be  entered   by  the  operator;    the  formula   used   is   similar  to  that 
developed   in   Subchapter  C  below.      The  computation  of  GDOP   is 
performed   in  the  receiver  by  the  following  formula: 


GDOP       =     [Trace  (A  A)-']'''  = 


(2.1) 


A  = 


cos 

ci 

COS 

4 

COS 

«1 

1 

COS 

4 

COS 

4 

COS 

"z 

1 

COS 

< 

COS 

4 

COS 

«^z 

1 

COS 

< 

COS 

"y 

COS 

< 

1 

COS   a^        =  direction   cosines  of  the  angles   between   the   range  to 
satellite   number   k   and  the   x   axis. 
An   acceptable   range  is   2  <  GDOP  <  4    (Jorgensen,    1984,    pp.    9-10). 


B.      PSEUDORANGE 

The  pseudorange  corrections  and  satellite's  clock  offset  are  calcu 
lated  and  applied  to  the  observed  pseudoranges  to  get  corrected  pseu 
deranges    (Milliken    &  Zoller,    1980,    pp.    4;    Spilker,    1980,    pp.    24). 


f\.bs=   f\ru*   c[^    ^    (5t^   -    5*3   )] 


obs 


tru 


c 


=   pseudo-range  observed   to   SV        [m] 

=  true   range   user  to  SV        [m] 

=  time  correction   consist  of  atmospheric, 

relativity,    and   earth    rotation        [s] 
=   SV     clock  offset  from  GPS  time        [s] 
=   user  clock  offset  from  GPS   time        [s] 
=  velocity  of   light        [m/s] 


(2.2) 


^\/yf^^ 


oJUJ^ 


17 


The  computation   of  satellite  clock  offset  from  GPS  time  is   done  by 
a   polynomial   model    (Equation   2.4),    given   by  Meyerhoff   (1985)    and  Van 
Dierendonck  et  al.    (1980).      The  model's   parameters   and  the  time   refer- 
ence of  the  parameters   are  given    in   the   navigation   messages.      The 
receiver's   clock  offset  from  the  GPS  time,    is   included   in   the  four 
unknown    receivers'   parameters. 

^tag~  *ts"   ^*sv  (2.3) 

6tsv             =  clock  offset  correction   of  SV  [s] 

^tag             ~   ^P^   time  tag   of  measurement  [s] 

tj-g               =  the  space  vehicle  code   phase  time  at  message 
transmission   time        [s] 


t 
SV  o        i      tag 


"-      =^o-   ^l^^agU^*   ^(W   toe)'  (2.4) 


a  =   parameter  of  clock   satellite  bias        [s] 

ai  =   parameter  of  clock   satellite  drift        [s/s] 

32  =   parameter  of  clock   satellite  aging        [s/s^] 

tgg  =   reference  time  of  clock  drift  at  GPS   system  time        [s] 

The   pseudorange  correction   consists   of  corrections   due  to   atmos- 
pheric,   relativity   and   earth    rotation   effects.      Formulas   for  these  are 
given    in   this   section.      The  atmospheric   correction    may   be   separated 
into  an    ionospheric  correction   and   a   tropospheric   correction    (Equation 
2.5).      The   ionospheric  correction   depends   on   the  electron   contents   of 
the  atmosphere.      The  best   ionospheric  correction   formula    uses   observa- 
tions  of   pseudorange   at  two  frequencies    LI    and    L2.      Equation   2.6   is 
used    in   the  computation.      The   ionospheric  correction   may   be  calculated 
by   an    ionospheric  model    included    in   the   satellite   navigation   messages; 
this   polynomial   model    is    less   accurate  than    Equation   2.6,    and   was    not 
used    in   this   evaluation.       If  the    P-code   is    not  available,    however. 
Equation   2.6  can't  be   used. 


^atm     '■ion*    '^  trop  (2.5) 
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The  ionospheric  correction  for  pseudoranges   is  given   by   Equation 
2.6   (Spilker,    1980,    pp.   25),    and  the  ionospheric  correction  for  pseudo- 
range   rates   by   Equation   2.7    (Meyerhoff,    1985,    pp.    20). 


tion=  1-5336  6t  g^u 


(2.6) 


T-on 
StgdLl 


=   ionospheric  correction        [s] 

=  difference  between   propagation   at   LI    and    L2    [s] 


*  ionrd 


(VQ^-./Q,) 


(2.7) 


[s] 


[(Ql/  Q2)'  -  ^   ]  ^ 

^onrd  ~   ionospheric   Doppler  correction 

Z^  =   LI    Doppler  count 

Q^  =   LI    frequency  multiplier  =   154 

y_  =   L2   Doppler  count 

Q2  =   L2  frequency  multiplier  =   120 

^  =   Nominal   satellite  frequency  =   10,230,000   Hz 

The  tropospheric  correction    is   computed   by  the  Chow  model 
tropospheric  correction   given    in    Equation   2.8    (Meyerhoff,    1985,    pp. 
22-24). 


t 


trop 


0.0276   P 


(2.8) 


[sin    e    (tan    9    +  0.0045)    +  0.00143] 

^  470  T^    E^-"   ^1.705  10^    p 

[T'    sin    e    (tan    0    *0.017)    *  0.00035] 

[mbar] 


E  =  35.65   H    10 

P  =   surface   pressure 

T  =  temperature        [°K] 

p  =  temperature   lapse   rate   set  at    .006        [deg/m] 

H  =   percent   relative   humidity 

9  =  elevation   angle  of   satellite        [deg] 

The  relativity  corrections    (tj-^-^)    due  to  the  eccentricity  of  the 
satellite's  orbit,    are  computed   by   Equation   2.9    (Meyerhoff,    1985,    pp 
27-28): 
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t       =     2  e  sin   E  y^^^   A^''^  (2.9) 

c 
t      ,  =   relativity  correction        [s] 

e  =  eccentricity  of  satellite  orbit 

E  =  eccentric  anomaly    (it  could   be  solved 

iteratively) 
y  =    WGS-72     value     of    the     earth's     universal     gravitational 

constant  =  3.986008  10^"        [mVs^] 
A  =   semi   major  axis   of  orbit        [m] 

c  =   speed  of   light   in   vacuum        [m] 

The  earth    rotation   corrections    (tj-^^-)    due  to  the  motion   of  the 
earth   during  the  signal's   propagation,    are  computed   by   Equation   2.10 
(Meyerhoff,    1985,    pp.    25-26]: 

trot=  ^^e   (   Y3^t>^     -   X33tYo)  (2.10) 

S^e  =  WGS-72  value  of  the  earth's    rotational    rate  = 

0.00007292115147        [rad/s] 
Xgat^at    ~   ECEF  of   satellite  coordinates,    WGS-72        [m] 
Xq  ,    Yq       =   ECEF  coordinates  of  the  station,    WGS-72        [m] 

C.      SATELLITE   POSITION   CALCULATION 

To  compute  a   ship's   position,    the   satellite   position   must   be 
computed   from   the  ephemerides   transmitted    in   the   navigation   messages. 
As   the   satellite  ephemerides   parameters   are   Keplerian   elements   in 
secular   drift  terms   and    harmonic  coefficients   only,    the   receiver  does    not 
need   sophisticated   software  to  compute  the   satellite  position.      It  must 
be   noted,    however,    that  the   parameters   change  every   hour.      The 
computed   mean   motion   was   solved   by    Kepler's   third   law.      The   Keplerian 
parameters    (mean   anomaly,    argument  of   latitude,    radius  of  orbital 
plane,    and   longitude  of  ascending    node)    are  affected   by   secular  drift 
terms   and   harmonic  coefficients,    which   must   be  corrected   from   their 
value  at   reference  time  by  their  drift.      The  drift  of  each    parameter   is   a 
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function  of  the  time  from  epoch  and   secular  drift  plus   harmonic  correc- 
tions.     A   sample  computation  of  satellite  position   may   be  found   in 
Mueller   (1964). 

1 .      Conversion  of  Receiver  Time  to  GPS  Time 


The  GPS  position   is  marked  by  the   receiver  clock;    however, 
the   reference  time  for  satellite  position   computation    is   GPS   system  time. 
The   receiver  time  tag   must  be  converted  to  the  GPS   system  time  by  the 
following   step. 

ttag=  W-    «V      -   ttran  (2.11) 

t^  =  GPS   system  time  at  time  of  transmission;    GPS   time  of 

tag  ^ 

reception   corrected   for  transit  time        [s] 

t  =   receiver  time  tag   at   receiver  clock        [s] 

rec  ^  •■    ■■ 

6t^  =   receiver  time  offset  from  GPS   system  time  originally   set 

to  zero,    it   is    included   in   the   unknown   parameter        [s] 
•ttran  ~  transit  time  of  signal   from   satellite  to   user        [s]  .  . 

The  time  from  epoch    is    related   to  the  ephemerides    reference  time  by: 

t       =  ttag   toe  (2.12) 

t  =  time  from  epoch,    if  greater  than  302,400   subtract 

604,800;    if  less  than   -302,400  add  604,800        [s] 
tog  =  ephemerides    reference  time  at  GPS   system  time    [s] 

2.      Computation  of  Orbital  Plane  Parameters 

The  orbital   plane  parameter   is   computed   by    Kepler's   third   law 
(Equation   2.13)    and    Kepler's   first   law    (Equation   2.16).      The  WGS-72 
parameters  for  the  computation   are: 

V  =  WGS-72   value  of  the  earth's   universal   gravitational 

constant  =  3.986008  10^"*        [mVs^] 
tl^  =  WGS-72  value  of  the  earth's    rotational    rate  = 

0.00007292115147        [rad/s] 
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The  mean   motion    is   affected   primarily   by   secular  drift  terms;    (Van 
Dierendonck,    1980,    pp.    67)    it   is   corrected   as  follows: 

n^     =     y^^^'  (2.13) 

n^  =  computed   mean   motion        [rad/s] 

A  =   semi   major  axis  of  orbit        [m] 

n       =   n     +   5n  (2.14) 

o 

n  =  corrected  mean   motion        [rad/s] 

5n  =  mean   motion   difference  from  computed   value 

M       =  Mq  -^   n   t  (2.15) 

M  =  mean   anomaly        [rad] 

Mq  =  mean   anomaly   at   reference  time        [rad] 

The  eccentric  anomaly  was   computed   by:  .  . 

M       =   E     *   e  sin    E  (2.16) 

E  =  eccentric  anomaly    (it  could   be  solved 

iteratively) 
e  =  eccentricity  of  satellite  orbit 

3.      Computation  of  Orbital  Plane  Parameters  in   ECEF 

The  orbital    plane   parameters    in    ECEF   are  affected   by   secular 
drift  terms   and    harmonic  coefficients;    these  effects   were  calculated   as 
corrections   to   the   parameters   as   follows: 

cos  V  =    (cos   E   -   e   )/(l    -   e  cos   E   )  (2.17) 

sin   V  =    [sin    E    (1    -   eM  ]    /    d    -   e  cos   E   )  (2.18) 

V  =  true  anomaly        [rad] 
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-    V 


(2.19) 


0  =  argument  of  latitude        [rad] 

u  =  argument  of  perigee        [rad] 

6U       =  C„^cos  20     *   C,^sin  20  (2.20) 

6U  =  argument  of   latitude  correction        [rad] 

Cyj,  =  amplitude  of  the  cosine   harmonic  correction   to  the 

argument  of   latitude        [rad] 
C  =  amplitude  of  the  sine   harmonic  correction   to  the  argu- 

ment of   latitude        [rad] 

6r       =   Cj-c  cos   20   +   Cj-gsin   20  (2.21) 

6r  =   radius   correction        [m] 

Crc  =  amplitude  of  the  cosine   harmonic  correction   term  the 

orbit   radius        [m] 
Cj-g  =  amplitude  of  the  sine   harmonic  correction   term  the  orbit 

radius        [m] 

61  =  Ci^cos  20   +   Ci^in  20  (2.22) 

51  =  correction   to   inclination        [rad] 

Cic  -   amplitude  of  the  cosine   harmonic  correction   term  the 

angle  of   inclination        [fad] 
Cig  =  amplitude  of  the   sine   harmonic  correction   term  the 

angle  of   inclination        [rad] 

U        =  0   +    6U  (2.23) 

U  =  corrected   argument  of   latitude        [rad] 

r       =   A    (1    -   e  cos   E)'^^    +    5r  (2.24) 


=  corrected  [m] 

=   L    ^   61  (2.25) 


o 


I  =  corrected   inclination        [rad] 

\q  =   inclination   angle  at   reference  time        [rad] 

Q       =  Q     *   ih     -  hj   t-  L  t  (2.26) 

o  o  e  e 

S2  =  corrected    right  ascension   of  ascending   node        [rad] 

fl^  =   right  ascension   at   reference  time        [rad] 

J^  =   rate  of  change  of   right  ascension        [rad/s] 

4.      Computation  of  Satellite  Position  in   ECEF 

Equations   2.27  and  2.28  are   used  to  compute  the  position   of 
the   satellite   in    its   orbital    plane    (Mueller,    1964,    pp.    147-170;    Van 
Dierendonck  et  al,    1980,    pp.    67): 

X^        =   r  cos   U  (2.27) 

Y'        =   r  sin    U  (2.28) 

X^  =   position    in   orbital    plane        [m] 

Y^  =   position    in   orbital    plane        [m] 

Equations  2.29   and   2.30  are   used   to  compute  the  position   of  the  satellite 
in    ECEF: 


X„,^=  X'   cos   9.   -   Y'   cos    I    sin   Q  (2.29) 


Y3^j.=  X'   sin   n   *   Y'   cos    I    cos   Q  (2.30) 


Zsat=   Y'    sin    I  (2.31) 


X        ,    Y        ,    Z  =   ECEF  of  satellite  coordinate,    WGS-72        [m] 

sat        sat  sat 
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D.      LEAST-SQUARES  ADJUSTMENT   (ITERATIVE  METHOD) 

The  corrected  pseudorange  between  an   unknown   station    P  and  the 
satellites   is   given   by  Jorgensen    (1984): 


Rcor   =   Rtru    *  c    (6t^   -    Stg^) 


(2.32) 


where: 


R 


cor 


R 


tru 


=  corrected   pseudorange  from   Satellite  Vehicle   (SV) 
number  k  to  the  unknown   station    P 

=  true  pseudorange  from   Satellite  Vehicle   (SV) 
number   k  to  the   unknown    station    P 


D  =  [(x,^,-x  y  ^  (Y33,-Y  y  ^  (Z3at-z  )'] 
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(2.33) 


D  =  the  preliminary  distance  between   SV  to  the 

unknown   station    P. 
XsatYsat^at=  the  coordinates  of  SV        [m] 

Xq,    Yq  ,    Zq  =  the  preliminary  coordinates   of  the   unknown    station 
P        [m] 
c  =  the  velocity  of   light   in   vacuum        [m/s] 

The  corrected   pseudorange  may   be   linerized   by   a   Taylor  series 
expansion   as: 


R^Qj.=    R^^-L  *6R|         AX  +  6RI         AY  +  6R|         AZ  +  6R| 


At 


6X|  6Y1  6Z|  6t| 

where  the  Jacobian   of   R   calculated   for   6X,    6Y,    5Z,    and    6t  are 


B   = 


6R 

6X" 


6R 

77 


6R 


6R 
"fit" 


Xo-Xs 


at 


or 


Y    -Y 


D 


o  'sat 

~D 
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where: 

R     ,  =  calculated   pseudorange  by  preliminary  coordinates  of 

P  from  SV  to  the   unknown   station    P        [m] 
The   least  squares  adjustment  may   be  done  by  the   iterative  method. 
The  observation   equation   matrix: 

AB    +  V  =  f  (2.34) 

f    =     R.ol     -     Kr^r-  (2.35) 

cal  cor 

where 

B  =  the  Jacobian   matrix,    n   by  4. 

A  =  the  correction   of  the   unknown   coordinates   matrix,    4  by 

1    matrix;    consists  of   6X,    AY,    AZ,    and   At 
V  =  the   residual   observation   matrix,    n   by   1 

f  =  the  difference  between    range  corrected  and    range  calcu- 

lated  by   a   priori   coordinates   matrix-,    n   by   1 
W  =  the  weight  matrix,    n   by   n 

w     =   1/a;      o   is   the  variance  of  the  observations 
n  =   number  of  observations,    minimum  4  observations 

R     -,      =  calculated   pseudorange  matrix   by   preliminary  coordinates 

of   P  from   SV  to  the   unknown   station    P 
f^cor      ~  corrected   pseudorange  matrix   SV  from  the   unknown 
station    P 
Minimizing   Equation   2.34  as   is   done  by  Mikhail   and   Ackerman 
(1974) 

N   =   B^WB  (2.36) 

t  =  B%f  (2.37) 

A=   N-^   t=   (b'^WB)"^    (B'^Wf)  (2.38) 

>^,    V,  ,    Zp     =  the  corrected   coordinates  of  the   unknown   station    P        [m] 
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Xp  =  >^  *  A(1)  (2.39) 

Yp  =  Yq  *  A(2)  (2.40) 

Zp  =  Zo  -^  A(3)  (2.41) 

6t^  =  A(4)  (2.42) 

Xq,    Yq  ,    Zq    =  the  preliminary  coordinates  of  the  unknown   station    P 
[m] 
A  =  estimator  of   reference  variance  matrix 

A(i)  =  element  of  matrix   A        [m] 

E.      KALMAN    FILTER 

Meyerhoff   (1985)    gave  a   method  for  the  computation   of  dynamic 
positions   by  means  of  an   eight-state   Kalman   filter  after  pseudoranges 
have  been   corrected   and   smoothing   has   been   applied  to  them   by   Doppler 
count.      The  computation   of  position   consists  of  computation   of  corrected 
pseudoranges,    satellite's   positions,    variance  of  pseudoranges,    Jacobian 
matrix    B,    Kalman   gain    (K),    updated   position,    and  the   updated   covari- 
ance  matrix    (P).      [For  details  of  the   Kalman   filter   see  Meyerhoff 
(1985),    Bierman    (1982),    and    Brown    (1983)].      The  variances  of  the 
pseudoranges   are  computed  from  the  variance  pseudoranges   of   LI    and 
L2.      The  variances  of  the  pseudoranges  of   LI    and    L2  are 

Vli     =   K2    -    Ki[      &L      n        ]    [   O-^   *   M  (2.43) 

lO'^Ll  lO^Ll 

Vli  ~  variance  of  pseudorange  data   for   LI    or   L2,    respectively 

at  timetag        [s^] 
Ki8-K2        =  variance  bias   factors        [s^] 
3  =  code   loop   noise  bandwidth 

n  =  number  of  trackers 

C^j^  =0.1    of  carrier  power  to   noise  density  of   LI 

C^j^  =0.1    of  carrier  power  to   noise  density  of   L2 
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Bj^  =  code  predetection   bandwidth 

The  variance  of  pseudorange  is   computed   by: 

V=   t^Ll*   (VLi*\t2)]/   m^/Q^V   -1]^  (2.44) 

Vpsd  ~  cariance  of  pseudorange        [s^] 

i;(l<,k)        =  V     ,    for  tracker   k,    respectively 
psd 

K       =   P(-)    B^CB   P(-)    B^    +   Z)-^  (2.45) 

K  =   Kalman   gain   matrix 

Z  =  variance  of  pseudorange  matrix 

B  =  Jacobian   matrix 

P(-)  =  The  covariance  matrix   before   updated 

of   this    data    point,    the    initial    value    is    the    initial    vari 
ance  of  X   ,    Y   ,    Z   ,    fit  ;   of  diagonal  are  zero 

P(*)    =   P(-)    -   K   B   P(-)  (2.46) 

P(+)  =  The   updated   covariance  matrix 

R^^o  =   R^Kc.-Floi  -c  t  ^„^  (2.47) 

omc  CDS      cal  cor 

^mc  ~   pseudorange   range   residuals        [m] 

R  Lg  =   pseudorange  observed        [m] 

R     -.  =   pseudorange  calculated        [m] 

t  =  correction   of  pseudorange        [m] 

cor  ^  »  L     J 

c  =   speed   of   light        [m/s] 

X(*)      =  X(-)    *    K(1)    R^^^  (2.48) 

Y(-)      =  Y(-)    *   K(2)    R^^^  (2.49) 

X(-)  =  X  coordinate  before   update        [m] 
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Z(*)     =  Z(-)    *   K(3)    Ronic  (2.50) 

6t(*)   =  6t(-)    *   K(4)    F\,^^/c  (2.51) 

Y(-)  =  Y  coordinate  before  update  [m] 

Z(-)  =  Z  coordinate  before  update  [m] 

6t(-)  =   receiver  clock  offset  before  update        [s] 

K(i)  =   row   i   of  matrix    K 

X(*)  =  The  updated  X  coordinate  [m] 

Y(  +  )  =  The   updated   Y  coordinate  [m] 

Z(  +  )  =  The   updated   Z  coordinate  [m] 

6t(*)  =  The   updated    receiver  clock  offset        [s] 
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III.    THE  TI-4100  TEST   IN   MISSISSIPPI    SOUND 


A.      FIELD  OPERATION 

The   U.S   Naval   Oceanographic  Office   (NAVOCEANO)    conducted   a 
test  on   Mississippi   Sound   in   May   1984  to  evaluate  the  performance  of 
the  TI-4100  on   a   marine  platform    (small   launch)    and   to  determine 
whether  the  TI-4100   is  adequate  for   hydrographic   surveying.      The 
launch   position   was   determined   by  a  GPS    receiver  and     four   Del   Norte 
Trisponders  on   shore  stations;    two  of  the  Trisponders  were   located 
above   second-order  conventional    North   American    Datum-1927    (NAD-27) 
bench   marks,    while  the  other  two  Trisponder  positions  were  determined 
by  means  of  the  TRANSIT   SYSTEM    (Dupont  and   Dunn,    1984,    pp.    11-49 
and    11-50).      One  TI-4100  was   on   the   launch   and  one  was   at  a   third 
NAD-27  bench   mark.      Three  first-order  bench   marks  were  observed   by 
a   TRANSIT    receiver  to  determine  the   local   datum   shift  from   NAD-27  to 
WGS-72    (Dunn,    1985). 

Data  for   1984  Julian   days    136  and   138  collected   during   the 
NAVOCEANO  test  were  arranged    in    standard   GPS   data   exchange  format 
(Scott  and    Peters,    1983,    pp.    3-23)    and   sent  to   MPS  for  study.      The 
GPS   shore  station   data   and   Trisponder  data  were   stored  on   a   single 
data   tape.      The  GPS   shore   station   data   were  acquired   every   12 
seconds,    and   the  Trisponder  data   were  acquired   every   second.      The 
launch    receiver's    raw  GPS   data    (consisting   of  pseudoranges, 
pseudorange-rates,    satellites   ephemerides;    timetag   and   atmospheric 
data)    were  acquired   and    recorded   every   second;    the   launch's    real-time 
positions   as   determined   by   GPS    (X,    Y,    Z  and    6t)    were  acquired   every 
6  seconds;    and   both   the   raw  data   and    real-time  positions  were  stored 
on   three  additional   tapes. 
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B.      DATA   PROCESSING 

NSWC  has  four  computer  programs  for  post-processing   raw  GPS 
data   into  the  receiver's  position   and  clock  offset.      The  first  program   is 
called  C0N9TR;   the  second   SOLTOM;   the  third   KALMN2;    and  the  fourth 
program   DOALL.      The  first  converts  data  from  the  standard  data 
exchange  format  to  a   new  format  compatible  with   the  second   and  third 
programs.      The  SOLTOM  and    KALMN2   programs   convert  the  output  data 
of  C0N9TR   into   ECEF  coordinates   by   least  square  adjustment  or  by 
eight-state   Kalman   filtering,    respectively.      DOALL  compares   dynamic 
post-processed  data  with   the  Trisponder  position   data   in    UTM  coordi- 
nates.     The  programs   were  written   by  Meyerhoff    (1985)    for  the 
Computer   Data   Corporation   Cyber-175  and   in   their  entirety  are  quite 
large.      I    was   not  able  to  convert  these  programs   for  the   NPS    IBM-3033 
due  to  their  complexity  and   size. 

The  data  from  GPS   time   185000  s   to   188600  s    (At  =  3600  s)   on    1984 
Julian   day   136,    was   converted  from   standard  9-track  tape  to  a    new 
format,    so  that   it  would   be  compatible  for  further  processing   by  the 
least  squares   adjustment  or   Kalman   filter  computer  program. 

For  the  3600-s   interval   for  which    shore  station   data   were  analyzed 
using  the   least  squares   adjustment  program,    the  positions  obtained   by 
post-processing   did    not  significantly   differ  from   real   time  positions. 

For  the  3600-s    interval   for  which   GPS   and   Trisponder   launch    posi- 
tions were  analyzed,    a   total   of  3171 -s   data   points   from   both   positioning 
systems  were  available.      The   launch   positions   determined   by  GPS 
consisted  of  2740-s   of  data   points   based  on   four  satellites   and  431 -s   of 
data   points   based  on   three   satellites.      The  positions  of  the   launch 
determined   by  GPS   were  post-processed    using  the   Kalman   filter  program 
and   compared  with   the   launch   positions   determined   by  Trisponders. 
Since  the  time  tags  of  the   shore  station's   GPS   positions   differ  from  the 
launch's,    further  differential   mode  post-processing  of  this   data   was    not 
done. 

The  post-processing   for  GPS   shore  station   positions   and   launch 
positions    used   broadcast  ephemerides.      Post-orbit  ephemerides   were   not 
available  to  possibly   improve  the  accuracy. 
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C.      ACCURACY  OF   LAUNCH    POSITION   DETERMINED  WITH 
TRISPONDERS 

The  "true"   position   of  the   launch   was   computed   by   least  squares 
adjustment  from  the  two  second-order  shore  stations    (NAD-27)    and  the 
two  TRANSIT   stations   in  WGS-72   using  the  four   ranges  from  the 
Trisponders.      The  two  second-order  shore  stations    (NAD-27)    were 
transformed  to  WGS-72  and  converted  to   UTM.      The  two  TRANSIT 
stations   in  WGS-72  were  converted  from   ECEF  to   UTM.      The  total   accu- 
racy of  the   launch's   position   determined   by  the  Trisponders   is   i   7  m,    a 
value  which   depends  on   the  accuracy  of  the  Trisponder  locations   and 
the  accuracy  of  the  Trisponders  themselves;    Table   I    shows   the  compu- 
tation. 


TABLE    I 

THE  TOTAL  ACCURACY  OF   THE   LAUNCH'S   POSITION 
DETERMINED   BY  TRISPONDERS 


The  accuracy  of  transformation  =  ±   6  m 

The  accuracy  of  Trisponder  =  ±  3  m      . 

The  accuracy  of  convertion   to   UTM  =  ±    1    m 

The  accuracy  of  Transit  =  ±    1    m 

The   rms   accuracy   of  the 

launch    position   aetermined 

by   Trisponder  =     ±   6.86  m 

The  accuracy  of  the   launch's   position   determined   by  the 
Trisponders   depends   on   the  accuracy  of  the   position   of  the  two   second- 
order   Class   2,    shore   stations    used   as   Trisponder   sites    locations.      For  a 
conventional   geodetic   station    Bossier   (1984)    gives   second-order  accuracy 
as   1:20,000.      The   length   of  a   side  is   generally   not   less   than   5000  m, 
and   the   standard   deviation   of  angle  observation   does    not  exceed   0.8". 
Thus,    the   position   accuracy  of  the   second-order  shore   station    sites    is 
approximately   ±   0.25  m.      The  accuracy  of  the  transformation   from 
NAD-27  to  WGS-72  was   determined   by   datum   shift.      The  datum   shifts 
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were  determined   locally   by  TRANSIT  observations  on   three  first  order 
NAD-27  stations.      The  root-mean-square  accuracy  of  the  transformation 
may  be  estimated  as  approximately  5  to  10  m   (Kumar,    1985).      The  accu- 
racy of  conversion  from  WGS-72  to  UTM  coordinates  may  be  assumed   ± 
1.0  m   (Kumar,    1985;    U.S  Army,    1973,    pp.    8-9) 

The  accuracy  of  the  launch's  position  determined  by  the 
Trisponders   depends   also  on   the  position   accuracy  of  the  other  two 
shore  stations  determined  by  TRANSIT  satellite  observations.      The 
accuracy  of  the  TRANSIT   stations   may   be  assumed   to  be   ±    1.0  m,    since 
observations  from  more  than   25  passes   were   used,    and   precise  ephemer- 
ides  were  applied    (Bossier,    1984,    pp.    3-10  and  3-11;    Kumar,    1985). 

Because  the  Trisponder   is   a   microwave  positioning   system,    it   is 
subject  to  multipath   errors   and    range   holes.      In   general   the  accuracy 
of  a   Trisponder  is   ±  3  m    (Umbach,    1981,    pp.    AD-2),    although    it  might 
be   larger  in    some  environments.      Because  the  Trisponder  observations 
were  made  at   night  during   satellite  passes,    the  accuracy  of  the  test 
data   may   be  better  than    ±  3   m. 

D.      ACCURACY  OF  THE   LAUNCH'S   POSITION    DETERMINED   BY  THE 
GPS   RECEIVER 

The  GPS   post-processing   position   values   were  compared   with   the 
launch's   position   determined   by  the  Trisponders,    and  the  mean   of  the 
rms   values   was  found   to   be   less   than    ±    1 1    m   if  four   satellites   were 
available   (Table   II).      Positions   for  the  first  three  seconds   were   poor, 
with   the  maximum   rms   difference   being  215  m.      The  mean   of  the   rms 
values  was  found  to  be   less   than    ±   21    m   if  only  three   satellites   were 
available   (Table   III). 

The  signal   quality  of   satellite   number  8  was   variable;    this   was 
indicated   by   bad   pseudorange   residuals,    which   caused   an    increasing   of 
the   northing,    easting   and    rms   differences.      The  signal   quality  of  satel- 
lite number  8  tended   to  be  cyclic,    decreasing    until    it  was   lost,    leaving 
only  3   satellites,    and  then    increasing   again,    making  four   satellites 
available  again.      The  signal   quality   is   measured   by  the  pseudorange 
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TABLE   il 

DIFFERENCE   BETWEEN   GPS-   AND  TRISPONDER-DETERMINED 
POSITIONS  WITH    FOUR   SATELLITES    IN    VIEW 


Mean 

Maximum 

Minimum 

Sigma 

rms 

10.677 

28.733 

0.422 

12.735 

Easting 

-7.058 

12.751 

-25.622 

10.653 

Northing 

2.136 

25.445 

-19.429 

6.854 

residuals,    which   are  considered   bad   when   exceeding   20  m.      A  total   of 
865  available  data   points   indicated   bad   pseudorange   residuals,    which 
were  the   result  of  bad   Doppler  counts.      The   number  of  bad   Doppler 
counts   were:      1    from   satellite   no.    6,    1    from   satellite   no.    3,    97  from 
satellite   no.    11,    and   766  from   satellite   no.    8. 

The  presence  of  the  four   "good"    satellites   and  four  good   pseudo- 
ranges  will   significantly   improve  the  accuracy. 

TABLE    III 

DIFFERENCE   BETWEEN   GPS-   AND   TRISPONDER-DETERMINED 
POSITIONS  WITH   THREE   SATELLITES    IN    VIEW 


• 

Mean 

Maximum 

Minimum 

Sigma 

rms 

20.538 

133.621 

0.765 

27.538 

Easting 

-15.920 

12.132 

-125.801 

25.193 

Northing 

-1.614 

32.296. 

-45.040 

11.120 

The  accuracy  of  the   launch's   position   determined  with   the 
Trisponders   is   better  than   the   launch's   position   determined  with   GPS. 
Although   29   percent  of  the  data    have   similar  accuracy,    further   studies 
to  separate  good   data   may   be  performed.      The  TI-4100  does    not, 
however,    have  the  distance   limitation   of  the  Trisponders:    It   is    limited 
only   by  the   present  coverage  of  GPS   satellites. 

The   real-time  dynamic   positions   were  compared   with    position   deter- 
mined  by   post-processing.      The  mean    position   differences   were  AX=  80 
m,    AY=   7   m,    and   AZ=   13  m. 
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IV.    CONCLUSIONS 

1.  The  overall   accuracy  of  Trisponder-determined  positions  for  this 
study  was   ±  7  m.      As  the  observed   range  increases,   the  accuracy  of 
positions  obtained  from  the  Trisponders  decreases.      Beyond   line-of- 
sight  Trisponders  cannot  be  used. 

2.  The  accuracy  of  ship   positions   determined   at   sea   by  a   TI-4100 
GPS   receiver  is  slightly  lower  than  the  position   accuracy  obtained  with 
shore-based   Trisponders,    although   29   percent  of  the  data   analyzed   have 
similar  accuracy.      The  GPS    receiver   is   more  versatile  than   the 
Trisponder  system   in   that  GPS   is   not   limited   by   line-of-sight.      Post 
processing   accuracy  may   be   improved   if  post-orbit  ephemerides   data 
were  available  for  the  analysis. 

3.  To   improve  the  position   accuracy   by   GPS    receivers,    further 
studies   in   a   differential   mode  should   be  performed   in   the  marine  envi- 
ronment. 

4.  The  main   difficulty  of  this   study  was   the  complexity  of  the 
programs  for  the  CDC  computer.      The  conversion   of  the   NSWC   computer 
program  to  the   MPS    IBM-3033   is   suggested  for  the  further  application   at 
NPS.      A  collection  of   hydrographic  subroutines   in   a   common    library  at 
NPS  would   be  very   helpful    in   future   hydrographic  studies   at   NPS. 
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APPENDIX  A 
ABBREVIATIONS  AND  ACRONYMS 


AFP  =  Approval'for   Full    Production 

C/A-code  =   Coarse  and   Acquisition   code 

CDC  =   Computer  Data   Corporation 

DSARC  =   Defense  System   Acquisition    Review 

Council 

DT&E  =   Development  Test   £•   Evaluation 

ECEF  =   Earth   Centered   Earth    Fixed 

FNOC  =    Fleet   Numerical   Oceanography   Center 

GDOP  =  GEOMETRIC   Dilution   of   Precision 

GPS  =  Global    Positioning   System 

IBM  =    International    Business   Machine 

lOT&E  =    Initial   Operational   Test  and   Evaluation 

JPO  =  Joint   Program  Office 
MHz                               .   =  Mega    Hertz 

NAD-27  =   North   American    Datum   1927 

NAVOCEANO  =   U.S   Naval   Oceanographic  Office 

NPS  =   Naval    Postgraduate  School 

NSWC  =   Naval    Support  Weapon    Center 

P-code  =   Precision   code 

rms  =   root  mean    square 

SEP  =   Spherical    Error   Probability 

SV  =   Satellite  Vehicle 

Timation  =  Time  and   Navigation 

USN  =   United   States   Navy 

UTM  =   Universal   Traverse  Mercator 

VLBI  =  Very-Long-Baseline- Interferometry 

WGS-72  =  World-Geodetic-System   1972 
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